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Abstract. We consider the problem of tidal disruption of a star by a super-massive rotating black hole. 
Using a numerically fast Lagrangian model of a tidally disrupted star developed in our previous works, we survey 
the parameter space of the problem and find regions where the total disruption of the star or a partial mass loss 
from the star takes place as a result of fly-by around the black hole. 

Our treatment is based on General Relativity, and we consider a range of black hole masses where the tidal 
disruption competes with the relativistic effect of direct capture of stars by the black hole. We model the star as 
a full polytrope with n = 1.5 with the solar mass and radius. We show that our results can also be used to obtain 
the amount of mass lost by stars with different stellar masses and radii. 

We find that the results can be conveniently represented on the plane of specific orbital angular momenta of the 
star (je,j<f,)- We calculate the contours of a given mass loss of the star on this plane, for a given black hole mass 
M, rotational parameter a and inclination of the trajectory of the star with respect to the black hole equatorial 
plane. In the following such contours are referred to as the tidal cross sections. 

It is shown that the tidal cross sections can be approximated as circles symmetric above the axis j$ = 0, and 
shifted with respect to the origin of the coordinates in the direction of negative je. The radii and shifts of these 
circles are obtained numerically for the black hole masses in the range 5 • 1O 5 M - 10 9 M o and different values of 
a. It is shown that when a — tidal disruption takes place for M < 5 • 1O 7 M0 and when a ~ 1 tidal disruption is 
possible for M < 1O 9 M . 



1. Introduction 

Tidal disruption of stars by a supermassive black hole 
may be very important in galactic centres. It may pro- 
vide enough gas to fuel central engines of AGNs, e.g. Hills 
(1975). It also may account for non-stationary flashes of 
radiation observed in certain non-active galaxies, see Recs 
(1988) for theoretical discussion and Komossa et al. (2004) 
for observations. 

Tidal disruption takes place when the periastron of the 
stellar orbit is smaller than the Roche tidal radius 

r T = {[^-R st « 1.5 • 10 13 f& tC m, (1) 
V m V m * 

where M is the black hole mass and to is the mass of star, 
R st is the stellar radius, M 7 = M/10 7 M Q , TO* = to/Mq 
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and Rsr = R s t/Ro,- Since the semi-major axis of the stel- 
lar orbit can be as large as a few pc, tidally disrupted 
stars have highly elongated (technically, parabolic) orbits. 
Another important process that competes with the pro- 
cess of the tidal disruption is the relativistic process of 
direct capture of stars by a black hole. The star can be 
directly captured by the black hole when its periastron is 
of the order of the black hole gravitational radius 

r g = M 3 . 10 i2 A f 7Cm . (2) 

Since for a black hole of mass larger than 10 7 M© the ratio 
r g/ r T > 0.2, the relativistic treatment of the process of the 
tidal disruption is important. 

We consider problem of tidal disruption of a star, or a 
partial stripping of mass from the star, by the black hole. 
Our treatment of the problem is fully based on General 
Relativity. We determine numerically an amount of mass 
lost by the star as a result of a tidal encounter, for differ- 
ent parameters of the stellar orbit and of the black hole. 
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Since the initial orbit is assumed to be parabolic and the 
initial state of the star is assumed to be unperturbed by 
tidal influence of the black hole, the terminology associ- 
ated with cross sections is used, see also Beloborodov et 
al. (1992) (BIIP) and below. As we have mentioned above 
(and will show quantitatively later on), the relativistic ef- 
fects are not important when the black hole mass is suf- 
ficiently small. On the other hand when the black hole 
mass is sufficiently large ~ 1O 8 M0 — 1O 9 M0 the stars can 
be directly swallowed by the black hole without tidal dis- 
ruption (see Section 3.1.1 and also BIIP and references 
therein). We obtain numerically approximate expressions 
for the amount of mass lost by the star for an intermediate 
mass range of the black hole, where the effects of General 
Relativity are important and the processes of tidal disrup- 
tion and tidal stripping of mass compete with the process 
of direct capture. 

The rotating black hole can be characterised by only 
two parameters, its mass M and the rotational parameter 
a = cLh/ {GM 2 ) < 1, where Lh is the black hole angular 
momentum. 

To describe the orbit of the star we must specify the 
orbit's integrals of motion. We introduce the spherical co- 
ordinate system (r,6,4>) with 9 = ir/2 at the black hole 
equatorial plane and consider the motion of the star very 
far from the black hole, r — ► oo. In this case a highly 
elongated orbit of the star can be characterised by the po- 
lar angle 9^ associated with the semi-major axis of the 
orbit, and by two projections of specific angular momen- 
tum, j0 and jg . By definition, these quantities are con- 
served during the motion of the star. Taking into account 
that the problem is symmetric with respect to reflection 
jif, — > — j0 we can consider only positive values of j^ > 0. 
Alternatively, one can characterise the stellar orbit by pro- 
jection of the specific angular momentum onto the axis 
of rotation of the black hole, L z and the square of the 
projection of the angular momentum onto the equatorial 
plane, Q. It is easy to show that the so-defined integral of 
motion Q coincides with the well-known Carter integral 
when the gravitational field of the black hole dominates 
over the gravitational field of the central stellar cluster. 
The quantities L z and Q are related to j^ and jg as 

L z = jo sin 600, Q = j\ + jg cos 2 6*00 . (3) 

The Carter integral Q is always positive for parabolic or- 
bits, and we will use its square root q = y/Q to have the 
same dimensions for the integrals of motion q and L Zl j^, 
jg. For 6>oo = 7r/2 and j^ > 0, L z = jg and q = j^. 

It is well known that for a given 9^ there is a region, 
S c in the plane (je,j<t>) corresponding to capture of the 
stars by the black hole (e.g. Chandrasekhar 1983, Young 
1976). By analogy, we define the cross section of the tidal 
disruption, St as a region in the plane (je,jcj>) where a 
partial tidal stripping of mass from the star or a full tidal 

1 Let the velocity components at the orbit apastron be 
and vg and the apastron distance be r a . Then, j$ = r a vg and 
je = r a v 4> . 



disruption of the star takes place (BIIP). An alternative 
definition of the tidal cross sections would be a region in 
the plane (L z , q) for a given 9^ where a tidal stripping or a 
tidal disruption takes place. We will call these as the cross 
sections in the plane (L z ,q). When 9 X — tt/2 the cross 
sections in the upper half plane {jg,jcf, > 0) and the cross 
sections in the plane (L z , q) coincide. It is very important 
to note that the cross section of direct capture of the stars 
in the plane (L z , q) does not depend on the angle 9^ (e.g. 
Chandrasekhar 1983) 2 . As we will sec (Section 3.1.3) this 
property is also approximately valid for the cross sections 
of the tidal disruption in the plane (L z ,q) (see the text 
below). This fact allows us to reduce the dimensionality 
of the parameter space and obtain the approximate cross 
sections of tidal disruption for 9^ ^ tt/2 by a geometric 
transform of the cross sections calculated for 9^ = ir/2. 

In order to calculate the cross sections of tidal disrup- 
tion one should evolve numerically a model of a tidally 
disrupted star in the relativistic tidal field of the black 
hole for a sufficiently long time. It should be done for 
different values j^, jg, 9^ as well as for different values 
of M and a. Approximately 10 3 — 10 5 runs are needed. 
Taking into account that the present day numerical 3D 
finite difference or SPH models arc rather time consum- 
ing it is impossible to survey the whole parameter space 
with these models 3 . Recently we have proposed another 
approximate model of a tidally disrupted star which is 
a one dimensional Lagrangian model from the point of 
view of numerical calculations, and therefore, it is numer- 
ically fast (Ivanov & Novikov 2001, hereafter IN, Ivanov, 
Chernyakova & Novikov 2003, hereafter ICN). This model 
is used in the present Paper. 

The model is briefly described in the next Section. 
Section 3 presents the results of our calculations. We dis- 
cuss our results and present our conclusions in Section 4. 

We use natural units expressing all quantities related 
to the star in terms of the characteristic stellar time t st = 
y i? 3 t /Gm, stellar radius R s t, etc. The parameters related 
to the stellar orbit are also made dimensionless. We use 
L z = L z /(cr g ), q = q/(cr g ), j^ = j^/{cr g ) and jg = 
je/(cr g ). We do not write tilde below assuming that all 
quantities mentioned above are dimensionless. 

We assume below that the star is a full polytrope with 
n = 1.5 and solar mass and radius. We show how to apply 
our results to the case of other values of the mass and 
radius. 

We mainly consider the tidal encounters of moderate 
strength, thus neglecting the possibility of violent tidal 
disruption and formation of strong shocks in the stellar 
gas. 



2 Technically, it follows from the fact that the radius of pe- 
riastron of a parabolic orbit does not depend on 9oo. 

3 See Ivanov & Novikov (2001) for an overview of works on 
tidal disruption and astrophysical applications. 
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2. The model of a tidally disrupted star 

The derivation of the basic dynamical equations of our 
model and detailed comparison of the results with results 
based on the 3D finite difference models, SPH model and 
affine models of the tidally disrupted star has been inten- 
sively discussed in our previous papers (see IN, ICN). 

We assume that summation is performed over all in- 
dices appearing in our expression more than once, but 
summation is not performed if indices are enclosed in 
brackets. The upper and lower indices describe the rows 
and columns of matrices. 

In our model, in order to reduce the number of de- 
grees of freedom, and accordingly, to significantly reduce 
the computational time we assume that the star consists 
of elliptical shells. The parameters and orientation of these 
shells change with time. This assumption allows us to ex- 
press the coordinates of some particular gas element in a 
frame co-moving with the centre of mass of the star, x l 
in terms of (non-evolving) coordinates of the same gas el- 
ement Xq in the unperturbed spherical state of the star 
(say, before the tidal field is 'switched on') as 



where the radius tq = \J xqiXq is the same for all gas 
element belonging to a particular shell, and e l = Xq/tq 
are direction cosines (eoie l = 1). We represent the position 
matrix and its inverse as a product of two rotational 
matrices A 1 , and E* , and a diagonal matrix : 



T) = A}B l m E™ 



V A]E\ 



(5) 



r(0 



and ai are the principal axes of the 



where B l m 
elliptical shell. 

The evolution equations for the matrix follow from 
the integral consequences of the exact hydrodynamical of 
motion: the so-called virial relations (e.g. Chandrasekhar, 
1969) written for a particular elliptical shell of interest, 
and have the form: 



T; = 3S?(9- l2nSi— {gSX 



-A k a k D k E + C k T , 

2 9 



Irrtk Tjln 



piny 



(G) 



where the dot stands for the time derivative. We use the 
mass enclosed in the shell corresponding to the radius tq, 
M = An J Q ° p§r\dr\ as a new radial variable, and po is 
the density of the star in its unperturbed spherical state. 
The matrix Sj represents the inverse to the position ma- 
trix Tj and g is the determinant of the position matrix. 
The matrix C\ represents the relativistic tidal tensor, and 
therefore it is symmetric and traceless. 

The dimensionlcss quantities D k have been described 
by e. g. Chandrasekhar (1969), p. 41. They have the form: 



du 



(T) 



where A = yj (a\ + u){a\ + u)(a§ + it). 

The quantities (£) and are determined by distri- 
bution of the density p and pressure p in the star: 



p Air 



dn-, 



and 

P tJ = ^ / 



(8) 



(9) 



where integration is performed over the unit sphere asso- 
ciated with the direction cosines e l . 

Obviously, the first two terms in equation [S] are due 
to action of the pressure force on the elliptical shell, and 
the last two terms describe the self-gravity of the star and 
action of the tidal forces. 

To complete the set of equations[H]wc should know the 
distributions of pressure and density over the star. The 
density distribution 



p(t,x l 



Po(ro) 
D 



(10) 



(4) 

is determined by the Jacobian D = | \ of the mapping 



between the coordinates x 1 and x 
as 



0; 



dx 3 

which can be written 



D{4) = 9 -^(S l m (Tn> + S h m (TD', 
z ' o 



(11) 



where prime stands for differentiation over rg. We use the 
standard relation p = kp 1 to find the distribution of pres- 
sure over the star. Here k is the entropy constant, and 7 
is the specific heat ratio. 7 = 5/3 is used later on. 

The dynamical equation |S] has the usual integrals of 
motion: the energy integral and the integral of angular 
momentum. Their form can be found in ICN. For our 
discussion it is important that in our model the energy 
integral can be naturally separated into kinetic, thermal 
and gravitational parts. In addition to these integrals, the 
quantities 



Xjk (M)=Titi 



(12) 



are conserved. They represent the conservation of circula- 
tion of the fluid over the elliptical shells in our model (IN, 
ICN). 

2.1. Numerical scheme and computational details 

Our numerical scheme is described in IN and ICN. The 
variant of the scheme used in the calculations is a non- 
conservative explicit numerical scheme, and we use con- 
servation of the integrals of motion to check the accuracy 
of the calculations. This simple scheme allows us to cal- 
culate a single tidal disruption event with a very small 
computational time. However, as was pointed out in ICN, 
this variant of the scheme suffers from a slowly growing 
numerical instability. A radical remedy would be an im- 
plicit conservative numerical scheme, but the schemes of 
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this type are much more time consuming. Therefore, in 
order to suppress this instability we average the dynami- 
cal variables over the neighboring grid points once every 
200 time steps. We have checked that this procedure leads 
only to a rather slow numerical leakage of the stellar en- 
ergy which seems to be non-significant for our purposes. 

We use a simple criterion of mass loss. We calculate 
the sum of kinetic and potential energy (per unit of mass) 
for each grid point and assume that when this sum is 
positive the corresponding elliptical shell is gravitationally 
unbound. This criterion is in good agreement with the 
results of 3D finite difference computations (see IN and 
ICN). 

In principal one can use other criteria of the amount 
of mass lost by the star. For example, one can consider 
the shells with positive total energy (i. e. the sum of the 
potential, kinetic and thermal energies) as being tidally 
stripped from the star. For a particular shell the ratio of 
the thermal energy E t h to the potential energy E p is pro- 
portional to a power of a characteristic size of the shell, 
R sh : E th /Ep oc R~h 3 ~ , ~ 4) - When 7 > 4/3 this ratio de- 
creases with time for an expanding shell. Therefore, tak- 
ing into account the thermal energy should not make a 
significant difference to our criterion. 

Our model cannot treat shocks which can, in prin- 
ciple, influence the criterion of the mass loss. However, 
we mainly consider tidal encounters of moderate strength 
where only a partial mass stripping of the star takes place 
and the presence of strong shocks is not expected. This 
statement is confirmed by results of 3D finite difference 
calculations where no evidence for a strong shock in the 
star has been found (e.g. Khokhlov et al. (1993), Diener 
et al. (1997)). Note that some shocks may develop after 
the tidal disruption event in a situation when partial mass 
stripping takes places. A gravitationally bound part of the 
stellar gas initially expelled from the star may return and 
hit the surface of a dense undisrupted stellar core pro- 
ducing a shock there. These shocks only redistribute the 
energy and momentum of the gravitationally bound part 
of the stellar gas and therefore cannot influence our crite- 
rion. 

For the calculation of the stellar orbits we use the usual 
Boyer-Lindquist coordinates (r, 9, <p) and the explicit form 
of the tidal tensor given in Diener et al. (1997) 4 . 

As follows from the results based on the finite differ- 
ence scheme (e.g Diener et al. 1997, and references therein) 
and our own results (IN, ICN) the stellar structure starts 
to evolve significantly only when the star approaches the 
periastron of the orbit, r p . Therefore we start our calcu- 
lations at the radius To = l-5r g and evolve our numerical 
scheme during the fly-by around the black hole until a fi- 
nal sufficiently large radius rfi n is reached. This radius is 
taken to be r fi„ = 4 — 5r p for a sufficiently small mass of 
the black hole, M < 8 • 10 7 A/q. In the opposite case we 

4 Note misprints in the Diener et al. (1997) expression for 
the components of the tidal tensor. The sign of all components 
should be opposite. 
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Fig. 1. The calculated cross sections of tidal disruption 
for the case a = 0.999 and different black hole masses. 
The levels of the amount of mass lost by the star, Mi ost , 
are shown by curves of different types corresponding to 
different masses of the black hole. Three curves of the same 
type show Mi ost = 0.1 (the outer curve), Mi ost = 0.5 (the 
middle curve) and Mi ost = 1 (the inner curve). Note that 
the curves corresponding to the different black hole masses 
may overlap (e.g. Mi ost (M = 0.1) = 1 and Mi ost (M = 
0.25) = 0.1). The solid curve shows the cross section of 
direct capture. Top: The case of small black hole masses 
is shown, M = 0.1, 0.25, 0.5 and 1. The triangle shows the 
result of 3D finite difference calculations of Diener et al. 
(1997). They obtained M% gt ~ 0.5. Bottom: The case of 
large black hole masses, M = 2, 4, 8, 16, 32, 64 and 100. 
Note that in the case M = 64 and M = 100 the curves 
corresponding to different Mi ost almost coincide, and we 
show only the curves corresponding to Mi ost = 1. 

terminate a single computation when two conditions are 
satisfied: 1) r fi n > 3r g ; 2) the dimensionless time r = i/i* 
from the beginning of computation is sufficiently large: 
t > 15. 

3. Results 

3.1. The case 600 = tt/2 

The bulk of our computations has been performed for 
6 '00 = 7r/2. In this case the tidal cross sections in the 
upper half plane (jg,j<p > 0) and the cross sections in the 
plane (L z , q) coincide. For symmetry reasons we can con- 
sider only the cross sections in the plane (L z ,q). As it is 
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Fig. 2. Same as Fig. 1 but a = 0.75 Top: M = 0.1, 0.25, Fig. 3. Same as Fig. 1 but a = 0.5 Top: M = 0.1, 0.25, 
0.5 and 1 Bottom: M = 2, 4, 8, 16, 32 0.5 and 1 Bottom: M = 2, 4, 8, 16 



shown later (see Section 3.3) the tidal cross sections in 
the plane (L z , q) almost do not depend on the value of 9oc 
and therefore the cross sections in the plane (j^,,je) can 
be obtained from those corresponding to 8^ = tt/2 by a 
geometric transform. 

3.1.1. The form of cross sections 

In order to calculate the tidal cross sections we specify the 
range of the black hole masses and rotational parameters. 
We calculate the amount of mass lost by the star, Mi ost , 
for a given value of the black hole mass and rotational 
parameter. These values are shown in Table 1. We plot 
the levels of constant Mi ost in the plane (L z , q) . The black 
hole mass M is expressed in units of 10 7 M Q and Mi ost is 
expressed in units of the stellar mass in all Figures. 

In Figures 1-4 we show the cross sections for o ^ 0. 
Three levels of Mi ost — 0.1, 0.5 and 1 are shown for a given 
a and M . Since the cross section of direct capture does not 
depend on M in our dimcnsionless units, we plot several 
different tidal cross sections and the cross section of di- 
rect capture corresponding to several different black hole 
masses in the same Figure. 

In the non-relativistic theory of tidal disruption the 
tidal cross sections are circles centred at the origin of the 
coordinates (L z ,q) with radius cx A/^ 1 / 3 in our units 5 . 



In physical units the characteristic size of the tidal cross 
section oc M 2 ^ 3 and the characteristic size of the capture cross 
section oc M. 



When the mass of the black hole exceeds a certain value 
the radius of the tidal cross section becomes smaller than 
the one corresponding to capture and the tidal disruption 
can no longer take place. In the fully relativistic treatment 
of the problem the situation is more complicated. We have 
checked with an accuracy of the order of 10 -3 that even in 
this case the levels of the mass loss always have a circular 
form. However, when a > the centres of these circles are 
shifted toward negative values of L z and the radii of these 
circles do not obey the simple oc M _1 / 3 law (see Figures 
7-11 and the Section 3.1.2). 

These effects are most prominent at a high rotational 
rate of the black hole. In Figure 1 we show the results of 
calculations for a = 0.9999. When the black hole mass is 
sufficiently small (see the upper part of this Figure) , the 
levels are approximately circular with a small shift toward 
negative values of L z . The triangle shows the position of 
the orbital parameters L z = 2.72945 and q = for the 
case intensively studied by Diener et al. (1997). They ob- 
tained Mi ost w 0.5 for M = 1.0853 x 1O 7 M . This value 
is in excellent agreement with our results for the mass 
M = 10 7 Mq. Note, however, that there is a considerable 
ambiguity in the results reported by Diener ct al. (1997) 
counting the gas elements leaving the computational box 
but having velocities less than the parabolic velocity as 
being stripped away from the star. Provided that these 
elements are not considered as being stripped, the value 
of Mi os t is considerably less, Mi os t ~ 0.32. The lower part 
of Figure 1 shows the case of high black hole masses. The 
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Table 1. The values of the black hole masses and rotational parameters used in the calcu- 
lations 

M/WMq O05 Ml 025 05 I 2 
M/1Q 7 M & 4 8 16 32 64 100 

a 0.25 05 0.75 0.9 0.999 



1 1 1 


1 , 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

M = 0.25 - 
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Fig. 4. Same as Fig. 1 but a = 
0.5 and 1 Bottom: M = 2, 4, 8 



0.25 Top: M = 0.1, 0.25, 



levels of the amount of mass lost by the star decrease to- 
ward the capture cross section with increasing mass. When 
M = M- ~ 4 • 10 7 M© the levels of mass loss intersect the 
capture cross section at the point (q = 0, L z = where 



(13) 



is a 'critical' negative angular momentum. It is defined 
for the equatorial orbits with negative angular momenta 
if the particles on the orbits with \L\ < |L_| are cap- 
tured by the black hole. When the black hole mass grows 
further, the intersection point moves toward the positive 
value of L z and when M > 1.6 • 10 8 M Q only the stars with 
positive angular momenta can be tidally disrupted. When 
1Q 9 M(t, the tidal cross section coincides with 



M = M. 



+ 



the capture cross section. In this case only the stars with 
angular momenta very close to 



L+ = 1 + Vl~ 



(14) 



can be tidally disrupted. Here the 'critical' angular mo- 
mentum L + is defined analogously to L_ for the or- 
bits with positive angular momenta (i. e. the stars with 
L < L+ are captured by the black hole). When M > M+ 
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Fig. 5. The amount of mass lost by the star calculated 
for the case a = 0. Since the cross sections in this case 
are concentric circles in the plane (L z ,q) we show Mi os t 
as a function of the radius L = W 'q 2 + L\ for the different 
black hole masses M = 0.1, 0.25, 0.5, 0.75, 1, 2 and 4. 
Vertical triangles show the position of the radius of the 
capture cross section L capt = 2. 




5 7 7.5 8 



Fig. 6. The amount of mass lost by the star compared 
with the results based on the Newtonian approach to the 
problem for the case a = O.The dash dotted curves corre- 
spond to extrapolation to the higher masses of the results 
obtained for M = 0.05 and the dotted curves take into ac- 
count the relativistic dependence of the periastron radius 
r p on L, see the text for details. 



only the process of capture of the stars is possible. The 
width of the region in the plane (L z ,q) where a partial 
stripping of mass takes place decreases with mass. As we 
will sec below this is mainly a relativistic effect. 

In Figures 2 — 4 we show similar results for smaller 
values of the rotational parameter, a = 0.75 (Figure 2), 
0.5 (Figure 3), 0.25 (Figure 4). In general, the behaviour 
of the levels of Mi ost looks similar to the previous case and 
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Table 2. The values of M_ and M + as functions of a 



a 


0.25 






0.5 






0.75 


0.999 


M_ 


~ 8 ■ 


1O Y M 


~ 8 


10 Y M Q 


< 8 


10 


M e 


- 4- 10 Y M Q 


M+ 


~ 8 ■ 


1O 7 M 


~ 1.6 


10 8 Mq 


~ 3.2 


10 ! 




~ 1O 9 M 



the levels of constant Mi ost can be well approximated as 
circles. The shift of the centres of these circles from the 
origin of the coordinates decreases with the decrease of a. 
The corresponding values of M_ and M + are shown in 
Table 2. As seen from this Table the value of M_ slightly 
increases and the value of M + decreases with decrease of 
a. 

When a = the cross sections are circles centred on 
the origin of the coordinate system, and we can use the 
dependence of M[ ost on L = L 2 . + q 2 to represent our 
results. This is shown in Figure 5. In this case M_ = 
M + ~ 4 • 10 7 M Q . Similary to the case of non-zero a the 
size of the region where a partial stripping takes place 
decreases with M. To show how important the relativistic 
effects are we compare our results with simple relations 
for the amount of mass lost by the star based on the non- 
relativistic theory of tidal disruption, see Figure 6. In the 
non-rclativistic theory Mi ost depends on the pcriastron 
distance r p , the black hole mass M, the stellar mass m* 
and the stellar radius R st only in combination 

V = (r P /r T f 2 , (15) 

where tt is the tidal radius given by equation ^ We 
calculate the dependence Mi os t(rf) for a very low mass 
M = 5 • 1O 5 M , and use this dependence to calculate 
Mi ost {L) for higher masses of the black hole (dot dashed 
curves in Figure 6). These curves give a good approxima- 
tion to the relativistic calculations for a small black hole 
mass M = 10 6 Af Q . However, the results diverge with in- 
creasing mass M . The curves based on the non-relativistic 
approximation systematically give a larger size of the re- 
gion where a partial stripping of the mass takes place. We 
have checked that when M > 10 7 M Q the deviation is very 
significant and these results are not shown in this Figure. 
To take into account some effects of General Relativity, 
Dicncr et al. (1997) have proposed using the relativistic 
dependence of r p on L in equation ED The curves calcu- 
lated according to this prescription are shown as dotted 
curves in Figure 6. The prescription of Diener et al. (1997) 
slightly improves the agreement, but the deviation is still 
rather large. This can be explained by the fact that there 
are relativistic effects not accounted for by the Diener et 
al. (1997) prescription, such as a different form of the tidal 
tensor, a different time spent by the star near the perias- 
tron, etc.. 

3.1.2. Semi-analytic representation of the results 
obtained 

As we discussed above the contours of the amount of mass 
lost by the star in the plane {L z , q) can be approximated as 




Fig. 7. The dependence of the shift of the approximating 
circles on mass M for three values of the rotational 
parameter a = 0.9999, 0.9 and 0.75. The lower and upper 
curves of the same type correspond to the circles approxi- 
mating the levels M[ ost = 1 and M[ ost =0.1, respectively. 

circles of a given radius tl shifted toward negative values 
of L z . This means that they can be represented in the 
form: 

(L z - df + q 2 = r 2 L , (16) 

where d < is the shift of the centres of the approximating 
circles. Therefore, we can describe the levels of Mi ost with 
two functions d(M, a) and r(M, a). We present these func- 
tions in Figures 7-11 for two values of Mi ost : Mi ost = 0.1 
and Mi os t = 1. Since in the non-relativistic approxima- 
tion the characteristic size of the cross sections scales as 
A/" 1 / 3 , it is convenient to introduce the rescaled radius 
and shift: d* = M 1 / 3 d and r* = M 1 / 3 Tl- The change of 
r* with mass is a purely relativistic effect. 

In Figure 7 we present the dependence of with mass 
for high values of the rotational parameter a = 0.9999, 0.9 
and 0.75. For any given curve the rightmost value of M is 
equal to M + . For a given value of a the shift corresponding 
to Mi ost = 1 is always larger than that corresponding 
to M^st = 0.1. In Figure 8 we show the same results 
calculated for the case of the small rotational parameters 
a = 0.75, a = 0.5 and a = 0.25. 
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M/10 7 M Q 

Fig. 8. Same as Figurc^but a = 0.75, 0.5 and 0.25. 




M/10 7 M Q 

Fig. 9. The dependence of the radius of the approximating 
circles r* on mass M for three values of the rotational pa- 
rameter a = 0.9999, 0.9 and 0.75. The curves correspond- 
ing to the same value of a and describing two different 
values of the amount of mass lost by the star Mi ost = 0.1 
and Mi os t = 1 intersect at M ~ 10. For black hole higher 
masses the radius of the circle corresponding to Mi ost = 1 
has a radius r* larger than for Mj ost = 0.1. 




10 20 30 40 50 60 



M/10 7 M a 

Fig. 10. The ratio / of the length of a segment of the ap- 
proximating circle lying outside the capture cross section 
to the total length of the approximating circle as a func- 
tion of the mass M. All curves correspond to Mi ost = 0.1 
and different values of a = 0.9999, 0.9, 0.75 and 0.5. 
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Fig. 11. Similar to Figure 10 but now the ratio /i of the 
length of an approximating circle lying outside the capture 
cross section to the total length of the unit of the tidal 
cross section and the capture cross section is shown as a 
function of the mass M. 
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In Figure 9 we show the dependence of on mass for 
a = 0.9999, 0.9 and 0.75. Note that all curves are very 
close to each other. The radius corresponding to the level 
Mi os t = 0.1 is larger than that corresponding to M; ost = 
1 only for sufficiently small masses M < M intersect ~ 
10 8 Af Q . When M = M int ersect the curves corresponding 
to the different levels intersect and when M > Mi nter sect 
the radius of the cross section corresponding to Mi ost = 1 
is larger than that corresponding to Mi os t =0.1. However, 
in the region outside the capture cross section the level 
corresponding to Mi ost = 1 is situated inside the level 
corresponding to M\ ost = 0.1 due to a larger negative 
shift. When M < Mi nter sect, the curves with different a 
coincide. When a < 0.5, we have M + < M intersect and 
the curves corresponding to smaller values of a coincide 
with the curves corresponding to the larger values of a. 
Therefore, we do not show these curves in Figure 9. 

In Figure 10 we show the ratio / of the length of the 
circles lying outside the capture cross section to the total 
length of the circles. The level Mi ost = 0.1 is shown for 
the cases a = 0.9999, 0.9, 0.75 and 0.5. For a particular 
curve, the mass of the black hole corresponding to / = 1 is 
approximately equal to M_ and the mass corresponding to 
the end of a particular curve is approximately M + . When 
a < 0.5, M- almost coincides with M + and the curves are 
almost vertical. They are not shown in this Figure. 

Figure 11 is similar to Figure 10, but here we show 
the ratio /i of length of the circles lying outside the cap- 
ture cross section to the length of the boundary of unit 
of the tidal cross section and the capture cross section. 
This quantity has a direct physical meaning. In a real as- 
trophysical setting where the super-massive black hole is 
embedded in a stellar cluster, the processes of tidal dis- 
ruption and capture of the stars typically occur in the so- 
called regime of the empty loss cone (e.g. Frank & Rees 
(1976), BIIP and references therein). In this regime only 
the stars with orbital parameters close to the boundary of 
the unit can be either disrupted or captured. So, the quan- 
tity /i gives the fraction of tidally disrupted stars among 
all stars with 6oo ~ tt/2 that have been destroyed by the 
black hole (i.e. either tidally disrupted or captured). 



3.1.3. The case 0oo ^ tt/2 

As follows from equation[3]the law of transformation from 
the coordinates (L Zl q) to the coordinates (je,j<j>) has the 
form 
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Fig. 12. The tidal cross sections calculated for ©oo = © + 
(shown by the circles) in comparison to the tidal cross 
sections calculated for Ooo = tt/2 for a = 0.9999 and 
Miost = 1- Note that the case Ooo = ©- is equivalent to 
the case shown, by symmetry. 
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Fig. 13. The tidal cross sections in the plane (je,j<f>)- The 
angle Ooo = tt/3 and a = 0.9999. 
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Fig. 14. Same as Figure 13 but Ooo = t/6. 



je = i^/cosOoo, 34> = yq 2 - L l cot 



(17) 



are allowed 6 , where 



Since the expression in the square root in equationEl mus t ©± = tt/2 ± arctan | — | 
be positive, only the values of Ooo m the range 



(19) 



0_ < Ooo < 0+ 



6 Note that these inequalities are also valid for the angle 
0(t) that changes during the fly-by around the black hole: 



(18) e_<e(r)<e+. 
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Fig. 15. Same as Figure 13 but ©oo = 7r/12. 
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Fig. 16. Same as Figure 13 but 6oo = 7r/24. 



In order to check the dependence of our results on ©oc 
we calculate several tidal cross sections in the plane (L z , q) 
for Goo = 0+ and compare them to those corresponding 
to Ooo = 7r/2. The results of comparison are presented 
in Figure 12. They show that this dependence is practi- 
cally absent. This may be explained as follows. When the 
periastron distance r p is sufficiently large, the effects de- 
termined by the rotation of the black hole, such as the 
dependence of the results on ©ooi are determined by high 
order corrections to the non-relativistic results. In the op- 
posite case of a star that has the orbital parameters very 
close to those corresponding to the capture, the orbit of 
the star near the periastron is close to a 'barrel-like' criti- 
cal orbit of constant r p . In this case, the angle O changes 
between two limiting values 0_ and 0+ near the perias- 
tron, and the dependence of the results on 6 M is averaged 
out (see also BIIP). Therefore, this dependence appears 
to be unimportant in the two limiting cases. The results 
shown in Figure 12 indicate that this dependence is also 
small in an intermediate case. 

Therefore, for a given value of ©oo, we obtain the tidal 
cross sections in the plane (je,j$) transforming the results 
obtained for O^ = tt/2 to the new coordinates according 
to the transformation law El The results are shown in 
Figures 13-16. One can see from these Figures that the 
cross sections becomes more symmetric with respect to 



the origin of the coordinate system with decreasing of ©oo- 
In the limiting case ©oo = the cross sections are circles. 

Let us consider the situation when the tidal cross sec- 
tion and the capture cross section intersect. The relative 
position of these cross sections is mainly characterised by 
their intersection point. Let L*"* be the value of the z 
component of the angular momentum corresponding to 
the point of intersection on the plane (L z , q) , and jg nt be 
the coordinate of the same point on the plane {je,j<t>)- 
The values of L™* and jg nt are connected by equation 17. 
As follows from this equation, when L™* < 0, j@ nt shifts 
leftward with decreasing of Ooo from 7r/2 to 0. In the op- 
posite case L* nt > 0, jg nt shifts to the right. As follows 
from Figure 1, we have L™ 1 w for M w 1.6 ■ 1O 8 M , 
and as we see from Figures 13-16 the point of intersection 
corresponding to this mass has jg « for all ©oo- 



4. Discussion and conclusions 

We have obtained the cross sections of the amount of mass 
lost by a star tidally disrupted by a rotating black hole in 
the range of black hole masses where the effects of General 
Relativity are important. Our results can be used in a sit- 
uation where the process of direct capture of the stars by 
the black hole and the process of tidal disruption compete 
with each other. With a high accuracy the tidal cross sec- 
tions in the plane of angular momenta (jg,j<f,) are circles 
shifted toward negative values of jg for n ^ 0. We have 
found the values of radius and the shift of these circles 
as functions of the black hole mass M for several values 
of the rotational parameter a. The ratio of characteristic 
size of the tidal cross section to the size of the capture 
cross sections decreases with the mass M. Since the shift 
of the capture cross section toward negative values of jg 
is always larger than the shift of the tidal cross sections, 
the cross sections intersect each other at a certain value 
of the mass. The intersection point moves toward positive 
values of jg and when the black hole mass is sufficiently 
large: M > M + {a), the process of tidal disruption is no 
longer possible. We estimate M + (a = 0) ~ 4 • 10 7 Mq and 
M+(a = 1) - 1O 9 M . 

We have used the numerical model of a tidally dis- 
rupted star developed in our previous work (IN, ICN). 
Approximately 5 • 10 3 tidal encounters have been consid- 
ered, the number at least an order of magnitude larger 
than all results obtained elsewhere. 

Our results would be helpful for detailed studies of 
tidal feeding of the central black holes by the stellar gas. 

To model the star we use the simple n — 1.5 polytrope. 
Since the stars are also often approximated by polytropes 
with larger values of n, it will be interesting to obtain 
the dependence of the results on n. Such a study lies be- 
yond the scope of the present Paper. We note, however, 
that the polytropes with larger n are more centrally con- 
densed, and therefore, it is more difficult to disrupt them. 
Accordingly, the mass M + (a) should be smaller for larger 
n. 



12 



P. B. Ivanov&M.A. Chernyakova: Relativistic cross sections of tidal disruption 



We consider the stars having solar mass Mq and solar 
radius Rq. However, our results can be used for stars with 
different masses m and radii R st - Indeed, the stellar mass 
and the radius enter in our calculations only in combina- 
tion t s t/t gr , where t st = \j and t gr = Therefore, 
when the stellar mass and radius are not the solar ones, 
the cross sections obtained for the mass M correspond to 
the re-scaled mass of the black hole 

M re = J^{^) m M. (20) 
V m Rq 

As we have shown there is a rather significant region 
in the plane {je,ji/>) where only a partial stripping of mass 
from the star takes place. The future of a part of the 
star that remains gravitationally bound after the tidal 
stripping depends on many factors. Let us suppose that 
the orbit of this remainder is unchanged after the tidal 
encounter. Taking into account that the central density 
of the remainder is smaller than the density of the un- 
perturbed star and the remainder has significant internal 
motions and rotation, one can suppose that the remain- 
der would be tidally disrupted during the next tidal en- 
counter. However, the processes of distant gravitational 
interactions with other stars of the central stellar cluster 
could increase the value of the orbital angular momentum 
of the remainder and it could be present in the cluster for a 
certain period of time. The observational detection of such 
remainders e. g. in the centre of our own Galaxy would be 
a very convincing test of the theory of tidal disruptions of 
stars by a black hole. 

Quantitative results obtained above should be taken 
with caution. From the results of comparison between our 
'old' variant of the model of a tidally disrupted star formu- 
lated and developed in IN and our 'new' advanced variant 
formulated and developed in ICN it follows that when the 
amount of mass lost by the star is large enough the 'old' 
model gives a larger value of the mass loss, and the differ- 
ence between these models is of the order of 30%. Also, the 
difference between the 3D hydro-simulations themselves 
can be as large as ~ 50% or even larger depending on 
particular numerical scheme and criterion for the mass 
loss (e.g Diener et al. 1997, and the discussion above). 
However, we believe that our qualitative conclusions are 
robust, and future work on improvement of the numerical 
schemes would mainly lead to corrections of the numerical 
values characterising the process of the tidal disruption, 
such as e. g. the value of M + (a). 

We neglect the possibility of formation of strong shocks 
in the star. Therefore, our results may be changed if future 
finite difference simulations will show that shocks play a 
significant role in the range of parameters we consider. 

We have calculated only one quantity related to the 
outcome of the tidal disruption event: the amount of mass 
lost by the star, Mi ost . As has been pointed out by Lacy et 
al. (1982) and Rees (1988) the dynamics of the gas lost by 
the star depends significantly on the distribution of this 
gas over the orbital energies. The gas elements with neg- 
ative orbital energies form an eccentric disc and can, in 



principle, accrete onto the black hole 7 . The gas elements 
with positive orbital energies may leave the gravitational 
field of the black hole. Therefore, in a future work we will 
be calculate the cross sections characterising the distribu- 
tion of the gas lost by the star over the orbital integrals of 
motion. Also, it would be interesting to calculate the cross 
sections for other quantities, such as the energy, internal 
angular momentum and average density of the gravita- 
tionally bound remainder. 
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